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We study the response of a superconductor with a strong spin-orbit coupling on an external 
magnetic field. The Ginzburg-Landau free energy functional is derived microscopically for a general 
crystal structure, both with and without an inversion center, and for an arbitrary symmetry of 
the superconducting order parameter. As a by-product, we obtain the general expressions for the 
intrinsic magnetic moment of the Cooper pairs. It is shown that the Ginzburg-Landau gradient 
energy in a superconductor lacking inversion symmetry has unusual structure. The general formalism 
is illustrated using as an example CePtaSi, which is the first known heavy-fermion superconductor 
without an inversion center. 
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I. INTRODUCTION 

Superconductors with unconventional, or anisotropic, 
pairing have remained one of the favourite and most- 
studied systems in condensed matter physics for more 
than two decades. Any superconducting material in 
which the symmetry of the pair wave function is different 
from an s-wave spin singlet, predicted by the Bardeen- 
Cooper-Schrieffer (BCS) theory, can be called "uncon- 
ventional" . From the initial discoveries of superconduc- 
tivity in the heavy-fermion compounds, the list of exam- 
ples has grown to include the high-T c cuprate supercon- 
ductors, ruthenates, magnetic superconductors, and pos- 
sibly organic materials. In contrast, such popular novel 
superconductor as MgB2, in which the order parameter is 
an s-wave singlet, is still "conventional" despite its many 
unusual properties. 

Although the pairing mechanism in most if not all un- 
conventional superconductors is subject to much debate, 
their behavior can be well understood using the symme- 
try approach, pioneered in Refs. 0, 0, 0J- I* 1 partic- 
ular, the intrinsic anisotropy and the multi-component 
nature of the order parameter lead to a variety of in- 
teresting magnetic properties, such as the internal mag- 
netism of the Cooper pairs, multiple phases in the vortex 
state, and the upper critical field anisotropy near T c not 
described by the effective mass tensor in the Ginzburg- 
Landau (GL) equations, for a review see, e.g., Refs. 0,O|. 

In most of the previous microscopic calculations of the 
magnetic properties of unconventional superconductors, 
the model of an isotropic band in a centrosymmetric crys- 
tal has been used. Historically, this has its origin in the 
fact that an unconventional Cooper pairing was first ex- 
tensively studied in the context of the superfluid 3 He, 
which is indeed an isotropic Fermi liquid with a weak 
spin-orbit (SO) coupling 6]. Although taking into ac- 
count a realistic Fermi surface anisotropy in a crystalline 
superconductor is not believed to cause any drastic qual- 
itative effects, it might lead to some considerable quan- 
titative changes compared to the parabolic band model. 
The SO coupling in crystals is usually taken care of by re- 
defining the basis of the single-electron states: instead of 



the usual Bloch spinors, the Cooper pairs are now formed 
by pseudospin eigenstates 0. Then the only significant 
change in the superconducting properties, compared to 
the case without SO coupling, is that the system is no 
longer invariant with respect to arbitrary SU(2) spin ro- 
tations, which alters the symmetry of the order param- 
eter in the pseudospin-triplet channel 0, • A detailed 
analysis of the temperature dependence of the upper crit- 
ical field, including the band anisotropy, impurity scat- 
tering, and sometimes the Fermi liquid corrections, has 
been done using the quasiclassical (Eilenberger) method, 
see e.g. Ref. |7fl and the references therein. A disadvan- 
tage of this approach is that it assumes a constant den- 
sity of states of the normal electrons near the Fermi sur- 
face and therefore fails to capture some contributions to 
the intrinsic magnetism of the Cooper pairs. Additional 
complications arise when a superconductor with SO in- 
teraction lacks an inversion center. In a nutshell, the 
symmetry analysis of superconducting phases should be 
modified if the SO coupling is strong, because the twofold 
degeneracy of the single-electron bands is now lifted al- 
most everywhere in the Brillouin zone, which makes it 
impossible to introduce pseudospin and also suppresses 
most of the pairing channels 0. 

Although most superconductors do have inversion 
symmetry, there are some exceptions. Early examples in- 
cluded such materials as VsSi [{j and H1V2 in which 
a possible loss of inversion symmetry is associated with a 
structural phase transition in the bulk of the crystal. The 
existence of superconductivity was later reported in fer- 
roelectric perovskites SrTi03 an d BaPb03-BaBi03 
[l2j |. It was pointed out in Ref. 01 that the surface 
superconductivity observed, e.g. in Na-doped WO3 [l4j| . 
is generically non-centrosymmetric simply because of the 
fact that the two sides of the surface layer are mani- 
festly non-equivalent. Possible effects of the absence of 
inversion symmetry in the layered high-T c cuprates were 
discussed in Refs. [l5j]. Very recently, superconductivity 
was found in non-centrosymmetric compounds CePtsSi 
[3 andUIr [13. 

This article is aimed to study the magnetic properties 
of a clean superconductor with arbitrary pairing sym- 
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metry and band structure, with or without an inversion 
center. We focus on the strong SO coupling limit, which 
is believed to be the case in many unconventional super- 
conductors, in particular the heavy- fermion compounds, 
because of the presence of elements with large atomic 
weights, such as U, Ce, etc. In contrast to the previ- 
ous works, the starting point of our calculations is an 
effective band Hamiltonian, which describes the dynam- 
ics of the Bloch electrons in a magnetic field 18]. The 
superconducting pairing is introduced using a BCS-type 
weak-coupling model, generalized for the case of an un- 
conventional pairing symmetry. We derive the GL free 
energy microscopically, which allows us not only to calcu- 
late the upper critical field, but also evaluate the intrinsic 
magnetic moment of the Cooper pairs in a crystalline su- 
perconductor. To the best of the author's knowledge, a 
microscopic derivation of the GL equations for a super- 
conductor lacking an inversion center, in the presence of 
an arbitrary SO coupling, has never been done before, so 
we fill this gap here. On the other hand, although some 
of our results in the centrosymmetric case are not new 
and can be found scattered in the literature, we found it 
instructive to treat both cases within the same general 
framework, which also highlights the important differ- 
ences between them. 

The article is organized as follows. In Section [HI we 
discuss the properties of the Bloch electrons in a mag- 
netic field in the normal state, and introduce the single- 
band effective Hamiltonian. In Section UTTl we study the 
properties of a strong SO coupling superconductor in a 
magnetic field near T c , derive the linearized GL equa- 
tions in the lowest order in B, and calculate the internal 
magnetic moment of the Cooper pairs, in both the cen- 
trosymmetric and non-centrosymmctric cases. In Section 
dVI we apply the general formalism to CePtaSi. Section 
El concludes with a discussion of our results. 



II. SINGLE-PARTICLE PROPERTIES 

To develop the necessary framework for the analysis 
of the superconducting properties, we first need to un- 
derstand how a uniform magnetic field affects the single- 
electron states in a normal crystal with SO coupling (with 
or without an inversion center). While for free electrons 
with a parabolic dispersion p 2 /2m the magnetic Hamil- 
tonian is obtained by simply replacing p with a gauge- 
invariant momentum operator p + (e/c)A (e is the ab- 
solute value of the electron charge), the case of band 
electrons should be treated more carefully. 

In zero field, the single-electron Hamiltonian has the 
form 



(1) 



where and c are the creation and annihilation opera- 
tors of band electrons with the wave vector fc, e l/ (k) is 
the quasiparticle dispersion in the z^th band, which takes 



into account all effects of the periodic lattice potential 
and the SO interaction, and J^ fc stands for the integra- 
tion over the first Brillouin zone. We assume that there 
is no disorder in the crystal, so that A; is a good quantum 
number in the absence of external fields. The Matsub- 
ara Green's function of electrons, defined in the standard 
fashion: 

G VlV2 {ki,T 1 ;k 2 ,T 2 ) = -{T T c kl ui(Ti)cl 2U2 (T2)}, (2) 

is diagonal with respect to both the band index and the 
wave vector: 



G v (k,ui r , 



I 



iuj n - e u (k) '' 



(3) 



where uj n = (2n + l)7rT is the fermionic Matsubara fre- 
quency (we use the units in which kg = 1). 

In the presence of a nonzero uniform magnetic field _B, 
Eq. is replaced by 



Ho = ^2 c LA( fc > B ) c kv, 



(4) 



where £ is the effective one-band Hamiltonian in the k- 
space 01- The main technical difficulty in the deriva- 
tion of Eq. |0J is that the corresponding vector potential 
A grows linearly as a function of r, leading to diver- 
gent matrix elements of the Hamiltonian with respect 
to the zero-field Bloch waves. As was first pointed out 
by Peierls these non-perturbative features can be 
taken into account by simply replacing the wave vector 
k in the zero-field band dispersion e L ,(k) by the gauge- 
invariant combination k + (e/hc)A(f), where f = iV& is 
the position operator in the fc-re pres entation. Later, this 
idea was elaborated in Refs. |2(J,|2l|], where it was shown 
that the Peierls Hamiltonian corresponds in fact to the 
zero-order term in the expansion of the general effective 
one-band Hamiltonian in powers of B: 



£ v {k,B) = e^+B^Ji^+B.B^K) 
where K is an operator in the fc-space: 

K = k + —A(r) = k + i— ( B x 

hc K ' 2hc V dk 



; (5) 



[here and below we use the symmetric gauge: A = 
{B x r)/2}. Since the components of K do not commute: 
[Ki,Kj] = —i(e/Hc)eijkBk, the order of application is 
important, so that £ is assumed to be a completely sym- 
metrized function of -JQ. This can be achieved, e.g., by 
representing the expansion coefficients in Eq. JSJ , which 
are periodic in k, in the form of a Fourier series over the 
lattice vectors R, and then replacing k — ► K to obtain 



iRK 



, etc. 



the operators e^(K) — ^2 R ^v{R)c 

If the electron bands are degenerate in zero field due 
to spin or pseudospin (see Sec. Ill Al below') . then the ef- 
fective Hamiltonian £ and all the expansion coefficients 
are 2x2 matrices. The Green's function corresponding 
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to the Hamiltonian Q is not diagonal with respect to fc, 
because the system is no longer invariant under lattice 
translations (it is still invariant though under the mag- 
netic translations which combine the lattice translations 
with gauge transformations). 

Although the explicit expressions for the expansion co- 
efficients in Eq. (JSj can be derived, at least in principle, 
using the procedure described in detail in Refs. [2l( , some 
important information can be obtained from general sym- 
metry considerations. The full symmetry group Q of the 
system in the normal state is given by a product of the 
space group and the gauge group J7(l). Assuming that 
there is no magnetic order in zero field and omitting the 
lattice translations , we can write Q = GxKxU(l), where 
G is the point group of the crystal, which may or may not 
include the inversion operation /, and K is the time re- 
versal operation. At non-zero B, the Hamiltonian (0J is 
invariant with respect to time reversal only if the sign of 
B (and of A) is also changed, which imposes the follow- 
ing constraint on the function £ : £ V (—B)K = £ U (B). 
In addition, the expansion coefficients must have certain 
transformation properties under the action of the point 
group elements, in particular, the band dispersion e(fc) 
must be invariant under all operations from G. 

Further steps depend crucially on whether or not there 
is an inversion center in the crystal lattice, which deter- 
mines the degeneracy of the zero-field bands. 



A. Crystals with inversion center 

If the crystal has an inversion center, then the bands 
are two-fold degenerate at each fc, because the Bloch 
states V'fc-i- = ^hv and ipk- = Klipkv have the same 
energy, belong to the same wave vector, and are orthog- 
onal. These states can be chosen to transform under 
the action of the space group operations similar to the 
spin eigenstates, in which case they are referred to as the 
pseudospin states Q. Thus the bands can be labelled 
by v — (n, a), where a — ± is the pseudospin projection. 
Focussing on a single band, we can omit the index n, and 
the effective band Hamiltonian (0 becomes 



£ a p(k,B) = e(K)8 a p - //,//, ,; K \a j>a 3 



(6) 



where aj are the Pauli matrices, and both the zero-field 
band dispersion e(fc) and the tensor fiij(k) are invari- 
ant under all point group operations. It is easy to see 
that this form of the effective Hamiltonian is compatible 
with all the symmetry requirements, in particular that £ 
should be Hermitian and K- and /-invariant. Indeed, the 
time reversal operator is K = (ia^Ko, where Kq is the 
operation of complex conjugation, which changes fc — * 
— fc. Therefore, we have [a 2 £(—k,—B)a 2 ]* = £{k,B). 
Also, £(—k,B) = £(k,B), because of inversion symme- 
try. In the limit of zero SO coupling, the usual Zeeman 
interaction term is recovered: /J-ij(k) — > [iB^ij^ where hb 
is the Bohr magneton. 



The Green's function © is a 2 x 2 matrix in the pseu- 
dospin space, which satisfies the following equation in the 
frequency representation: 

(iu n - £i) Q7 G 7/3 (fci,fc 2 ;w„) = 6 Q/ 3<5(fci - fc 2 ). (7) 

The Fourier transform of the Green's function, defined 
as 



G a p{r u r 2 ;uj n ) = ^ e lkl 

fei ,k 2 



ri — ik 2 r 2 



G Q(3 (fci,fc 2 ;w„), (8) 



satisfies the equation 

(iuj n - £ 1 ) ai G 1 f3{r 1 ,r 2 \uJn) = 5 a p8{r\ - r 2 ), (9) 

where £ is the Fourier transform of the effective band 
Hamiltonian ©, which is obtained by simply replacing 
K by the real space operator 

^=-4 + ^ (r)= -4 + 2^ (Bxr) - 

The subscript 1 in £\ or £\ means that the operator 
acts on the first argument of G. It should be noted 
that the Green's function JSJl is not the same as the 
Green's function of the band electrons in the coordi- 
nate representation. The latter is defined as {r\o\(iio n — 
Ha)~ 1 \r 2 cj') = {ria\k 1 a)G al3 {k ll k 2 ;uj n ){k 2 l3\r 2 a') (the 
summation over repeated indices is implied), where 
(ra\ka) = tpkaifcr) is the Bloch spinor, with a =t,| 
being the z-projection of spin. 

The second term in K presents some difficulty because 
it grows linearly as a function of r and therefore cannot 
be treated as a small perturbation. To handle this prob- 
lem, we seek solution of Eq. in a factorized form 



G a p(r u r 2l uj n ) = G Q/3 (ri - r 2 , uj n )e 



i¥j(ri,7- 2 ) 



(10) 



where <p(r\, r 2 ) = (e/hc) A(r)dr, and the integration 
goes along a straight line connecting r\ and r 2 [z| • Using 
the identities 

[ 2 A(r)dr = TA(r h2 ) + hBx(r 1 -r 2 )} ) (11) 
c?n,2 J ri 2 

one can show that the translationally-invariant function 
G(ri — r 2 ) — G(R) obeys the equation 



{iuo n — £) ai Gjp(R,uj n ) = 5 a pS(R), 



(12) 



where the operator £ is obtained by replacing K in 
the argument of £ in Eq. @ by Kr — —id/dR + 
(e/2hc)(B x R). 

The advantage of introducing the function G is that, in 
contrast to Eq. @, the magnetic field term in Eq. 112|) 
can be treated as a perturbation at small enough B. The 
precise condition can be easily obtained in the case of an 
isotropic parabolic band e(fc) = h 2 (k 2 — k^)/2m, when 
the solution of Eq. i|12|) hi zer0 ne ld is G a /3(R,uJn) ~ 
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5 a pe 



where vp = hkp/m is the 



Fermi velocity. Because of the fast oscillations of the ex- 
ponential, the characteristic scale of the derivative d/dR 
is kp. On the other hand, the scale of R is given by 
Hvp/ksT, so that the field-dependent term in KrG is 
small compared with the gradient term if huj c -C fc^T, 
where lo c = eB / mc is the cyclotron frequency. Although 
this condition does not have a simple form for a realistic 
band structure, it is usually assumed that the perturba- 
tive treatment of B in Eq. I|12|l is legitimate at all but 
very low temperatures, where the Landau level quantiza- 
tion effects become important. 

The Fourier transform of G satisfies the equation 



£(k,B)] G 7/5 (fc,c 



(13) 



which is solved perturbatively in B. The expansion of 
the effective band Hamiltonian has the form 



£ a p(k, B) = e(k)8 aP - Bm af3 (k) + Q{B 2 ), (14) 



where 



m l:a p(k) = i — 



v(k) x A 



Sa/3 + ^i]{k)Oj,al3- (15) 



The first term comes from the expansion of e(K), with 
v(k) = (l/h)de(k)/dk being the band velocity, while the 
second one is obtained by replacing K with k. As ob- 
vious from Eq. (|14|) . m can be interpreted as the mag- 
netic moment operator of the band electrons, although 
one cannot say that the first and the second terms cor- 
respond to the orbital and the spin magnetic moments 
respectively, because both v(k) and fJ-ij{k) include the 
effects of SO coupling. The solution of Eq. (|13fl can be 
written as G = Go — BGomGo + 0(B 2 ). Inserting the 
expression (|15fl here and keeping only the corrections of 
the first order in B, we have 



G a f3(k,UJ n ) 



Sat 



iu n - e(k) 



-BiHijik)- 



iu> n ~ e(fc)] 2 



(16) 



Note that because of inversion symmetry, G Q ^(— fe, ui n ) 
G a p(k,u> n ). 



B. Crystals without inversion center 

In the absence of inversion center in the crystal lat- 
tice, the electron bands are non-degenerate almost ev- 
erywhere, except from some high-symmetry lines in the 
Brillouin zone. The formal reason for this is that without 
the inversion operation /, one cannot in general construct 
two orthogonal degenerate Bloch states at the same k 
(note that the Kramers theorem still holds: there is a de- 
generacy between the time reversed states -0fc„ and Kip^ 
belonging to k and —k respectively). The above is not 
valid at zero SO coupling. In that case, there is an ad- 
ditional symmetry in the system - the invariance with 



respect to arbitrary spin rotations, which leads to the 
bands being two-fold degenerate because of spin, so that 
the results of the previous section apply. 

Assuming that the SO coupling is strong and the bands 
are well split (which is the case in CePtaSi |23j), the 
effective single-band Hamiltonian Q can be written in 
the following form 



£{k,B) = e{K) - B\{K) 



(17) 



where the band dispersion e(fe) is invariant with respect 
to all point group operations, and A(fc) is a pseudovec- 
tor, which, being a property of the crystal in zero field, 
satisfies the conditions (g\)(g~ 1 k) = \(k), where g is 
any operation from the point group. Because of the 
time-reversal symmetry, we also have e(— k) = e(k) and 
A(— k) = — X(k). At a nonzero B we have £(~k, —B) = 
£(k,B), but £(—k,B) ^ £(k,B) in general, because of 
the lack of inversion symmetry. An example of the micro- 
scopic calculation of A(fc) using a simple two-dimensional 
model is given at the end of this subsection. Also, in Sec- 
tion llVl below. we discuss how to find the momentum de- 
pendence of A in a non-centrosymmetric tetragonal crys- 
tal. 

The only modification to the analysis of Sec. Ill Al is 
that both the effective Hamiltonian JSJ and the Green's 
function (J2J become scalar functions. The Green's func- 
tion is factorized: 

G( ri ,r 2 ;u> n ) = G{ ri -r 2 ,uj n y^' r -\ (18) 

where the Fourier transform of G satisfies the equation 

[iw n -£{k,B)]G{k,u n ) = l. (19) 

As in the centrosymmetric case, at low fields we solve 
this equation perturbatively in B, using 



£ (fe, B) = e(k) - Bm{k) + Q{B 2 ), 



where 



m(k) = i 



2c 



v(k) x 



Ok 



A(fe) 



(20) 



(21) 



has the meaning of the magnetic moment operator of the 
band electrons. The contribution from the first term in 
m to G vanishes, and we finally have, in the first order 
in B, 



G(k,u n ) = 



1 



BiXi{k), 



1 



\lUJ n 



e{k)f 



(22) 



Most of the previous works on non-cent rosy mmetric 
superconductors, both two-dimensional [13, Ea, H i an d 
three-dimensional J25|, have been based on the Rashba 
model (we would like to mention, in particular, Ref. 
[26^. in which the GL functional was derived for a one- 
component s-wave order parameter in a Rashba super- 
conductor). In this model, the combined effect of the SO 
coupling and the lack of inversion symmetry is mimicked 
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by an additional term in the single-particle Hamiltonian: 



H = '^to(k)a\ iy a ka +^'^2n-((j aa iKk)(J kry a k< j>. (23) 



Here c, a' =T,I is the z-axis spin projection, the opera- 
tor Ofco- destroys an electron in a Bloch state of energy 
eo(fe) corresponding to zero SO coupling, and n is a unit 
vector allowed by symmetry (in a 2D system, n is sim- 
ply the normal vector to the plane). Choosing n = z, 
we diagonalize the Hamiltonian (|23|l by a unitary trans- 
formation afccr = Uk,anCkn (n — 1,2), which gives two 
Rashba bands: 



ei(2)(fe) = eo(fc) ± |7l*i 



(fcj_ = \ k% + k%), with the eigenfunctions 



^fc,l(2)( r ) = 



= JLf 



y/2 V T» 



(24) 



(25) 



where tan^/- = k y /k x . The bands (|24|l are non- 
degenerate almost everywhere, touching only at the two 
poles of the Fermi surface along the z axis. We would 
like to emphasize that the band indices n = 1,2 cannot 
be interpreted as the pseudospin projections. Indeed, un- 
der time reversal the pseudospin eigenstates would trans- 
form similar to the spin eigenstates, i.e. into one an- 
other. However, being a symmetry of the Hamiltonian 
time reversal transforms the Rashba bands into them- 
selves, which can be directly verified for the eigenstates 
(2J : 



ie 



V2 



V2 
1 



Mi 



'■'■Pk 



-ikr 



e lkr oc V-fc,i, 



and similarly for ipk,2 (we used (f- k = <fik + t)- 

It is easy to show that in the presence of a non-zero 
magnetic field the effective Hamiltonian for the Rashba 
model can be cast in the form l|17|) . To obtain the pseu- 
dovector A(fe), let us consider a two-dimensional system 
in a field parallel to the xy plane. Then the Hamiltonian 
(|23[) is modified by the Zeeman term: Hb = Hq — /ib^B. 
The diagonalization of Hb, followed by an expansion in 
powers of B, gives 

£ 1{2) (k,B) = eo(fc) 



± \/l 2 kl + 2 lt i B (k x B) x + fi 2 B B 2 
e 1{2) (k) - X m (k)B, 



where 



^l(2)(k) = ±^B 



k x n 



(26) 



In this article, we want to keep our discussion as general 
as possible and therefore do not resort to any explicit 
model, such as the Rashba model, to describe the SO 
coupling. Our results are based only on the symmetry 
considerations and valid for an arbitrary strength of the 
SO coupling and any band structure. 



III. MAGNETIC RESPONSE IN THE 
SUPERCONDUCTING STATE 

A. Crystals with inversion center 

Now let us take into account the attractive interaction 
between the band electrons in the Cooper channel. The 
total Hamiltonian is given by H = H a + H int , where the 
free electron Hamiltonian H is given by Eq. Q and, for 
a BCS-type mechanism of pairing, the interaction part 
can be written as 



Hint = I Vr «/3, 7 5(fc,fe / )4+g/2,a c i 



fc+q/2,/3 



k,k' ,q 



XC -fc'+q/2,7 C fc'+q/2,<5- (27) 



The pairing potential does not depend on the external 
magnetic field and is assumed to have a factorized form: 



1 d 

V a0 , 7S (k,k') = -■=Vj2**,«(>W*l„s(. k ')' 



(28) 



o=l 



with the coupling constant V > 0. Here x P a (fc) are the 
2x2 matrix basis functions of an irreducible representa- 
tion r of dimensionality d of the symmetry group of the 
system at zero magnetic field . The pairing interaction 
is nonzero only inside a thin shell of width e c (the cutoff 
energy) in the vicinity of the Fermi surface e(fc) = 0, i.e. 
^a(k) — 5' a (fei?)/ c [e(fc)], where kp is a wave vector at 
the Fermi surface and the cutoff function / c (e) is local- 
ized about the origin, e.g. f c (e) = 9(e c — \e\). The basis 
functions are assumed to be orthonormal: 

^<tr[*t( fc )* 6 (fc)]) ( 

= \ <tr[*t(fc)* 6 (fc)]) / c 2 (c) = <W c 2 (e), (29) 

where the angular brackets denote the averaging over the 
constant energy surface e(fc) = e: 



((•••)>. 



1 



N (e] 



(30) 



and iVo(e) = J2k ~ e ( k )] ^ s the normal-metal density 
of states (DoS) per one pseudospin projection. 

It follows from anti-commutation of the fermionic op- 
erators that ^ a ,/3a(— k) = —^a^apik). In the presence of 
inversion symmetry, the even in k (pseudospin-singlet) 
and odd in k (pseudospin-triplet) pairing states can be 
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considered separately. In the singlet case, the matrix ba- 
sis functions can be represented in the form 

^a,af3(k) = {ia 2 )a(3^a{k), (31) 

where <f> a (k) are the even scalar basis functions of the T 
representation. In the triplet case, we have 

*a : Q/j(fc) = (ma 2 ) a p(pa,i{ k ) ( 32 ) 

where <p a (k) are the odd vector basis functions of the T 
representation 

The superconducting order parameter can be repre- 
sented as a linear combination of the basis functions: 

(k), (33) 

a 

with the coefficients r/ a playing the role of the order pa- 
rameter components, which determine, for instance, the 

I 



free energy T of the superconductor. In the vicinity of 
the critical temperature T C (B), one can keep only the 
quadratic terms in the expansion of T: 

F = Y J jdrr,l(r)S ab r }b (r). (34) 

ab 

Here S is a d x d matrix differential operator of infinite 
order: 

S ab = - fdR S ab (R)e~ mD , (35) 



where D = — iV r + (2e/Hc)A, and the translationally- 
invariant function S ab {R) is expressed in terms of the 
Green's functions (JI35J. Its Fourier transform is given by 



S ab (q) = T^T^A 

n k 



%is{ k ) G Pi (k + |, u/„) G a s y-k 



Q 

r 



(36) 



^a'f3js(k) 



= — exp 



4hc 



d d 
dki dk 2 



*L/3( fc l)*b,75(fc2) 

k\ —ho —h 

5*^(fc)*i,7«(*) + *^- fi (V fc < Q/3 x V fe * 6 , 75 ) + Q(B 2 ) 



8hc 

r 



The derivation of these formulas is outlined in Appendix 
TKl As obvious from Eq. I|35|l . the operator S is a com- 
pletely symmetrized function of the components of _D, 
which do not commute: [Di,Dj\ = —i(2e/hc)eijkBk- 
Also, its Taylor expansion contains only even powers of 
D, because S ab (—R) — S ab (R) due to the inversion sym- 
metry. 

The field dependence of the phase transition tempera- 
ture at arbitrary B can be found from Eq. (|34(l : T C (B) 
is defined as the temperature at which the minimum 
eigenvalue of the operator S passes through zero. For 
an isotropic s-wave order parameter, the corresponding 
equations were derived and solved in Refs. |27l . w hile 
for an isotropic p-wave case it was done in Ref . |28j . In 
a general case, i.e. for an arbitrary band structure and 
pairing symmetry, T C (B) can only be calculated numer- 
ically. 

Here we focus on the properties of our superconductor 
in the weak field limit. We have T = j F dr, where the 
free energy density can be represented as 

F = A abV * a r, a + A„,, , ,,,/.;/),/),,///, - MB. (37) 

This expression has the usual form expected on the phe- 
nomenological grounds, with K ab ^j being the general- 



ized effective mass tensor, and M having the meaning of 
the intrinsic magnetic moment of the Cooper pairs. The 
linearized GL equations follow from Eq. (|37|l after the 
minimization over the order parameter: SJ-/5rj2{r) = 0. 
Below we outline how to calculate the free energy density 
using our weak-coupling model. 

The first term in F is obtained by putting q = B = 
in Eqs. (|35|l and (|36() . which gives 



A 1 r 

A ab — y°ah 



i^tr[*t( fe) ^ (fe)]4S[e(fe)]j 



(38) 



and 



1 , e 
— tanh — . 

2e 2T 



The necessary momentum cutoff in Eq. (|38|) is provided 
by the basis functions ^ a (fc), which are restricted to the 
vicinity of the Fermi surface. Calculating the momentum 
integral with the help of the normalization condition l|29|l . 
we obtain A ab = [(1/V) — I]S ab , where 



I(T) = / de JV (e)£(e)S(e) ~ N F ln 



2e c 



irT 



(39) 
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(C ~ 0.577 is Euler's constant). To obtain this re- 
sult we made the usual assumption that N (e) is a 
slowly-varying function within the energy shell of width 
e c near the Fermi surface, which allows us to replace 
it by a constant - the DoS at the Fermi level Np — 
Nq(0). At the zero-field critical temperature T C) we 
have I{T C ) — 1/V, which gives the standard BCS re- 
sult: T c ~ 1.13e c exp(— 1/NpV). Expanding A ab in the 
vicinity of T c , we recover the familiar expression for the 
uniform term in the free energy density: 



A ab = a(T - T c )S a b 



(40) 



where a — Np/T c . 

Next, we calculate the intrinsic magnetic moment 
M. Using the small- B expansions of the normal-state 
Green's function G and the vertex A, we obtain in the 
singlet case: 



M, 



Afro 



where ((■■•))o stands for the Fermi-surface averaging (j20J, 
and I is defined by Eq. To derive this expression, 

we again used the fact that the basis functions are non- 
zero only in a narrow vicinity of the Fermi surface, which 
allows one to separate the energy integration from the 
integration over the Fermi surface. A similar calculation 
in the triplet case gives 

Mi = ^(MxVA),) I 



+2*77*77;, ^y(fe) (0* x 4> b )^h, 



where 



h(T) = J de N {e)f c {e)S 1 {e) ~ — ^ In (41) 
5 1 (e)=T^-J- 

^ ' I If, __ — i 



1 _ 1 dSje) 
-iu) n — e 2 de 



Here N' F = Aq(0) is a measure of the electron- hole asym- 
metry near the Fermi surface. Putting T — T c , using the 
BCS result for the critical temperature, and choosing real 
basis functions (which can always be done if the normal 
state is non-magnetic) we finally obtain the density of 
the intrinsic magnetic moment of the Cooper pairs: 



where j ab = -~/ ba is given by 



in the singlet case, and 



f d<j> a {k) 8^ b {k) \ 

4HcV Bijl \ dkj dh / 



(42) 



(43) 



e 1 



d(j) a , m (k) d(f>b, m {k) \ 
dki / 



' N F V 



1 " e jk l {llij(k)<j>a t k(k)<pb,l( k ))o ( 44 ) 



in the triplet case. It follows from these expressions that 
M = for any order parameter corresponding to a one- 
dimensional representation of the point group, both in 
the singlet and triplet cases. 

Finally, let us evaluate the gradient terms in Eq. 137f) . 
The magnetic field dependence of the coefficients K ab .ij 
can be neglected, which follows from the fact that the 
lowest eigenvalue of the operator K ab ^jDiDj is already 
linear in |_B|, see Appendix iBl The physical meaning of 
this is simple: the suppression of the critical temperature 
due to the gradient energy is always linear in a weak field, 
regardless of the dimensionality of the order parameter 
and the shape of the Fermi surface. Taking the second 
order derivative in Eq. (|36|l at B — and calculating the 
Matsubara sums, we obtain: 

K ab ,ij = -\& (tT[^ a (k)^ b (k)]v i (k)v j (k)) I 2 
-h 2 (tr[*t( fe )^( fe )] m -i( fe )) o/l . 

Here m~ 1 {k) — (l/h 2 )d 2 e(k)/dkjdkj is the inverse ten- 
sor of effective masses, I\ is defined by Eq. (|41|) . and 



h(T) = J deN (e)fZ(e)S 2 (e) 
where 



7C(3) 
8tt 2 T 2 



N f , (45) 



1 



1 



1 

2 {iuj r , 

1 

16T 2 e 



(iuj n 
1 



ty -iui n - e 
1 



-ILOr, 



sinh ( -— - ) cosh 



27V 



\2TJ 



has a peak near e = 0, and £(s) is Riemann's zeta- 
function. 

Putting all the pieces together, replacing T with T c , 
and using real basis functions, we finally have 

K ab>ij = 7 -^§!^N F (Mk)Mk)v l (k)v 3 (k)) 



8V N F 
in the singlet case, and 
7((3)h 2 



h 2 N' 

+ (Uk)Mk)m^(k)) Q (46) 



16tt 2 T, 2 



Nf ((f) at i(k)(f) b! i(k)vi(k)vj(k)). 



'8VN F 



<^, ; (fc)0 M (fe)m- 1 (fe)) o (47) 



in the triplet case. Assuming a spherical Fermi sur- 
face, a completely isotropic pairing corresponding to 
the unity representation of G, and neglecting the 
electron-hole asymmetry, Eq. (|46|) yields Kij = 
(5 ii [7C(3)n 2 /487T 2 T 2 ]A F u 2 , j2|. For an anisotropic Fermi 
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surface, but still a conventional pairing, the results of 
Ref. |23 are recovered. 

Now we would like to make a few comments about 
our results. The internal magnetism of superconductors 
has been discussed before mostly for a charged isotropic 
Fermi liquid without SO coupling, see, e.g. Ref. j^D]. In 
this case, the density of the pair magnetic moment can be 
divided into the orbital and spin parts, both being small 
due to the smallness of both the quasiclassical param- 
eter (kp^o) 2 *C 1 (£0 is the coherence length), and the 
electron- hole asymmetry N' F Q. Here we do not make 
any assumptions about the strength of the SO coupling, 
therefore the orbital and the spin magnetic moments can- 
not be separated, in general. For a general band dis- 
persion, one can neglect neither of these contributions 
apriori, before calculating the Fermi-surface averages in 
Eqs. i|43|l and i|44|) . In particular, the energy dependence 
of the single-electron DoS in the metals with d- and /- 
electrons is usually quite significant, which can lead to 
an appreciable electron-hole asymmetry near the Fermi 
level. 

In terms of the response of the superconductor on a 
weak external field, the gradient terms produce a linear 
in B suppression of T c , see Appendix [BJ The value of 
the slope dH&jdT can be calculated either analytically 
(in very few cases), or using a variational approach. On 
the other hand, the pair magnetism can compete with 
the gradient energy, leading even to the possibility of 
increasing T c as a function of B, if the internal magnetic 
moment is large enough. Such mechanism was recently 
proposed in Ref. [3lJ to explain the phase diagram of the 
ferromagnetic superconductor ZrZn2. 

B. Crystals without inversion center 

In this case, the calculations are somewhat simpler be- 
cause the bands are non-degenerate. We assume that 
the Cooper pairing occurs only between the electrons in 
the states with opposite momenta, which are transformed 
into each other by time reversal. Then the most general 
BCS-type Hamiltonian can be written in the form 

Him = H^t + H int + H inl ( 48 ) 

where 

H int = lET, V n 1 Hk,kyt n ct kn C- k , nCh , n 
n k,k' 

Hint ~ 2 Ktm^jk ) C kn C -kn C -k>mCk'm 

n^m k.k' 

Hint — 2 ^imft ^') c fe„ C -fc.m C -fc'«i C fc'»- 

n^m k.k' 

Here n and m label the non-degenerate single-electron 
bands, e.g. the Rashba bands l|24|) . The Hamiltonian 

H^ t describes the intra-band pairing, describes the 



pair scattering between the bands, which can result in the 
superconducting gaps induced on more than one sheet of 

(3) 

the Fermi surface, and H int corresponds to the pairing 
of electrons from different bands. 

A considerable simplification occurs if the supercon- 
ducting gaps are much smaller than the interband en- 
ergies. For example, the band structure calculations of 
Ref. [2^| show that the SO band splitting in CePtaSi ex- 
ceeds the superconducting gap by orders of magnitude. 
In this situation, the formation of interband pairs de- 

(3) 

scribed by H\ n [ is strongly suppressed for the same rea- 
sons as for the paramagnetically limited singlet supercon- 
ductors [32] : the interband splitting cuts off the logarith- 
mic singularity in the Cooper channel, thus reducing the 
critical temperature. Although the bands may touch at 
some isolated points at the Fermi surface, as is the case 
for the Rashba bands (|24l) at k || z, the interband pairing 
in the vicinity of those points is still suppressed due to 
the phase space limitations. We also neglect the possi- 
bility of the Cooper pairs having a non-zero momentum 
(Larkin-Ovchinnikov-Fulde-Ferrell phase) ($3, which is 
expected to be suppressed as well by the large depairmg 
effect of the SO band splitting. 

In this paper, we further neglect the inter-band pair 

scattering process described by H in ' t , leaving the investi- 
gation of its effects for future work. Thus, we focus on 
a single non-degenerate band for which the pairing be- 
tween time-reversed states |fe) and K\k) ~ | — k) near 
the Fermi surface can be written as 

Hint = 2^2 ^( fe ' k ') c k c Kk c Kk'C k ', (49) 

k,k' 

where <^ Kk denotes the creation operator of an electron 
in the state K\k), and the pairing potential is assumed 
to have a factorized form 

d 

V(k,k') = -VY / MkWa(k')- (50) 

a=l 

with V > 0. Here (f) a (k) are the scalar basis functions 
of an irreducible representation T of the point group of 
the crystal in the absence of magnetic field, which are 
nonzero only inside the energy shell of width e c near the 
Fermi surface: 4> a (k) = </> a (fcj?)/ c [e(fc)], and orthonor- 
mal: 

(#(fc)&(k)) e = mk)Mk)) Q / c 2 (e) = <W c 2 (e). (51) 

The parity of the basis functions can be determined us- 
ing the following arguments |3J|. Although the time- 
reversed state K\k) belongs to the wave vector — k, it 
is not the same as | — k). In fact, K\k) = t(k)\ — fc), 
where t(fe) is a non-trivial phase factor, which satisfies 
t(— fe) = — t(fe). This allows us to write c^ Kk = t(k)c_ k 
and CKk = t*{k)c—k- Inserting these relations in Eq. 
(|4*9"j) . we have 

Hint = \ V(k, k')c{cl k c^k>ck> , (52) 
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where V(k,k') = t[k)t*(k')V(k,k'). From the anti- 
commutation of fermionic operators it follows that 
V(k, k') has to be an even function of both arguments, 
i.e. one should choose even basis functions <fi a (k) in 
the expansion l|50[) . Treating the interaction (|52(l in the 
mean-field approximation, one obtains the order param- 
eter A(fe) = t{k) J2 a Va<Pa{k), which is odd in k. In Ref. 
|23j . the nodal structure of A(fe) was analyzed in terms 
of the odd basis functions. This has been corrected in 
Ref. [341, where the the importance of the phase factor 
t(k) was recognized. 

Allowing for the possibility of a non-uniform supercon- 
ducting order parameter, the Hamiltonian (|52|l becomes 



Hint — - V (^' fe ') C l+q/2 C -fc+q/2 



k,k F ,q 



The order parameter can be represented as 

A(k,q) = J2va(q)*a(k), (54) 



where ^aik) = t(k)<f> a {k) = —^ a (—k) satisfy the or- 
thonormality condition (4 r *(fe)4 r b(fe)) c = 5 a bfc( e )> see 
Eq. (E3J. 

The contribution to the free energy quadratic in the 
order parameter has the form (|34(l with the kernel now 
given by 



(55) 



xc_ 



k'+q/2 c k'+q/2- 



(53) where S a b(R) is the Fourier transform of 



Sab(q) 



T^5> Qb (fc)G(fc+^ n jG 

n k 



2' 



2 



(56) 



Aafc(fc) 



exp 



Ahc 



d 



d 



dki dk-2 



*S(fei)*6(fca: 



*:(fe)* 6 (fe) + 1— B(V k V* a x V fc tf 6 ) + 0{B 2 ). 
4nc 



r 



The derivation is similar to the centrosymmetric case, see 
Appendix El 

An important difference from the previous case is that, 
although the functions x f' a (fc) still have a definite parity, 
the Green's functions (|2"2"|l do not: G(— fc, u> n ) ^ G(k, w n ) 
in general, therefore S ab (—R) 7^ S ab (R)- This means 
that the expansion of the free energy density now con- 
tains gradient terms of an odd degree in D: 



F = C vlm + fZ vtDiVb + fZi V* a D i D jV b + ..., (57) 
where 

tib = ^ ab ~ S ab (q = 0), 



f(2) 



f (l) 
J ab.i 



f (2) 
J ab 7 ij 



etc. 



dS ab (q) 



q=0 



1 d 2 S ab (q) 

2 dqidqj 



q=0 



Using Eq. (J2H, it is easy to see that = at B = 0. 

Keeping only the lowest order terms in the free energy 
density expansion in a weak field, we have: 



F = A ab r]*rjb + K a b,ij r{ f a D l Djr\b - MB 
+K a b,ijBi ri*Djr]b, 



(58) 



where K a b,ij = d ff£j / 1 dBi\ T =T a ,B=a- The uniform con- 
tribution to F can be calculated in the same fashion as 
in the previous section, and we obtain 



A ab = a(T — T c )5 ab , 



(59) 



where the critical temperature T c is given by the same 
BCS expression as in the centrosymmetric case, but now 
a = A F /2T C . 

The pair magnetic moment M and the generalized 
effective mass tensor K ab ^ can be calculated similarly 
to the centrosymmetric case. Using real basis functions 
4> a {k), we obtain 



M 



e 1 
8&U' 



and 



7cm 2 



32ir 2 T? 



N F ((f, a (k)<f>b(k)v i (k)v j (k)) Q 



(60) 



h 2 N' F 
16V Np 



(Mk)Mk)m^(k)) . (61) 



To calculate the coefficient K ab ,ij, we expand S ab (q) to 
the first order in both B and q and evaluate the Mat- 
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TABLE I: The character table and the examples of the basis 
functions of the irreducible representations of C^. 



which 



r 


E 


Caz 


fi 


even 4>r{k) 


odd cf>r(k) 


A 1 


1 


1 


1 


k x ~\~ ky ~\~ ck z 


k z 


A 2 


1 


1 


-1 


(kx ky^k x ky 


(kx ky}k x kyk z 


Bi 


1 


-1 


1 


m x n,y 


{kx ~~ ky)k z 


B 2 


1 


-1 


-1 


kx ky 


kx kyk z 


E 


2 








kxkz j kyk z 


kx j ky 



subara sums, which gives 



K ah>ij = h(r a (k)Mk)^(k)v J (k)) I 2 



where A(fe) is the momentum-dependent pseudovector 
that determines the linear response of the band electrons 
on a weak magnetic field, see Eq. (|17ll . and lip are de- 
fined by Eqs. (|41(l and 1|45|) respectively. Using real basis 
functions, we finally have 



E~ab,ij 



7 ^^N F (Mk)Mk)\(k)v 3 (k)) Q 



1 



N' F 



<f) a {k)(j) b (k) 



d\j(k) \ 
dkj /„' 



(62) 



Note that the phase factors t(k) have dropped out of both 
K a b,ij and K a b^j. To evaluate the Fermi-surface averages 
in Eqs. H60I61I62|) explicitly, one has to know the band 
structure [including A(fc) and t(k)] and the momentum 
dependence of the order parameter. 



IV. APPLICATIONS TO CePt 3 Si 

CcPtsSi is a heavy- fermion material without inversion 
center, which was recently found to become supercon- 
ducting at T ~ 0.75K jjjjj. It has a tetragonal lattice 
symmetry described by the point group G = C^, which 
is generated by the rotations C^ z about the z axis by 
an angle 7r/2 and the reflections a x in the vertical plane 
(100). The Fermi surface is invariant under all the op- 
erations from dv and also the inversion, the latter be- 
ing the consequence of the time-reversal symmetry. The 
band structure calculations of Ref. [2^ show that the SO 
coupling in this material is strong and therefore the de- 
generacy of the bands is lifted everywhere, except along 
the z axis. 

The point group has five irreducible representa- 
tions: four one-dimensional (A\, A 2l B±, and B 2 ), and 
one two-dimensional (EYsee Tabled Although the order 
parameter is odd in k |23| . its nodal structure is deter- 
mined by the even basis functions |34j | . Here we consider 
only the case of a one-component order parameter, for 



A(fc,r) = J7(r)*(fc) = r)(r)t(k)<f>(k), 



(63) 



where <p{k) = <p(—k). The pair magnetic moment van- 
ishes, and the GL free energy (|58|l takes the form 

F = a(T — T C )M 2 + Kij^DiDjT) + A , , />',//' /),;/. 

Dropping the terms proportional to N' F and using the 
symmetry of the Fermi surface, we have 



K., = = ^ = i^S^F <^(fc)^(fc)) , 



K, 



32tt 2 T 2 



327T 2 T2 ' 



N F {4?(k)vl{k)) . 



(64) 



In order to calculate Kij, we need an expression for 
A(fe), which satisfies the following symmetry require- 
ments: A(-fe) = -A(fc), (C^AXC^/fc) = A(fc), and 
(a x \)(a~ 1 k) = A(fc) (since A is a pseudovector, we have 
ct x \ = IC 2x X = C 2x X, where C 2x is a rotation by an 
angle 7r about the x axis). To solve these constraints, we 
represent A as an expansion over the odd basis functions 
of the irreducible representations of , see Table [I] 



A ( fc ) = EE A rAa(fe), 



(65) 



r 0=1 



where (j){—k) = —<j)(k). It is straightforward to check 
that only the representations A 2 and E contribute to the 
expansion i|65|) . so that the most general expression for 
A(fe), which satisfies all the symmetry requirements, is 
given by 

A(fe) = X E <t>E,2(k)x - 4>E,i(k)y + \A 2 <pA 2 (k)z, (66) 

where and Xa 2 are constants. Substituting it into 
Eq. (|62|l . using the fact that the Fermi velocity v(k) 
transforms according to a vector representation E + A\ , 
and dropping the terms proportional to N' F , we finally 
have 



E^ X y E\y X K 



Km , 

8n 2 T? ' 



N F U\k)<i>E,x{k)v x {k) 



(67) 



All other vanish by symmetry. 

Finally, the GL free energy density can be written as 



F = a(T-T c )\r,\ 



■rf [K 1 {Dl + Dl)+K 2 Dl\r } 

+K V *{B x D y - B y D x ) V . (68) 



While the second-order gradient terms here are typical 
for a one-component order parameter in a uniaxial crys- 
tal, the last, linear in both D and B, term is unusual and 
occurs only because of the absence of inversion symmetry. 

As an application of the above results, let us calculate 
the upper critical fields for B parallel and perpendicular 
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to the z axis. To this end, we solve the linearized GL 
equation obtained from Eq. I|68|l . If B — B(0, 0, 1), then 



(69) 



If B = £>(cos tp, sin^, 0), we choose the gauge A = 
Bz(sin cp, — cos (p, 0). The lowest eigenvalue of the GL op- 
erator corresponds to the order parameter with no mod- 
ulation along the field direction: 



r](r) = exp 



i— (B x r) z z 
he 



where zq is an arbitrary parameter. The function f(z) 
satisfies an equation which can be reduced to the stan- 
dard harmonic oscillator equation by a shift in the coor- 
dinate: z = Z + zq + (hcK / AeKi). Thus we find 



Ginzburg-Landau theory is applicable. Starting from the 
effective single-band Hamiltonian in the magnetic field, 
we obtained the expressions for the GL effective masses 
and the internal magnetic moments of the Cooper pairs 
in terms of the Fermi-surface averages, for an arbitrary 
pairing symmetry and crystal structure, both in the cen- 
trosymmetric and non-centrosymmetric cases. 

For a superconductor without inversion symmetry, un- 
usual terms, linear in both the magnetic field and the 
order parameter gradients, were found in the free energy 
expansion. The order parameter itself corresponds to the 
pairing of electrons in the time-reversed states within the 
same non-degenerate band. As a simple application of 
our general formalism, we derived the GL functional for 
CePtsSi. It was found that although the unusual gradi- 
ent term does not affect the upper critical field in a bulk 
sample, it could result in a field-induced enhancement of 
T r in a thin film. 



and the field-dependent critical temperature 



T C (B) = T C (B = 0) 



2e VK1K2 
he a 



B 



K 2 



4aK 



(70) 



B 2 , (71) 



which is completely isotropic in the xy-plane. We see 
that, surprisingly, the K-term does not affect the linear in 
B suppression of T c , giving rise only to a small, quadratic 
in field, correction to T C (B). Neglecting the latter effect, 
we find 



H c 2,xy{T) — 



he 



:(T C -T). 



(72) 



The last term in Eq. (|71ll could become dominant in a 
film of CePtaSi. If the thickness of the film is less than 
the correlation length £ 2 = K2/a(T c — T), then the order 
parameter (|70|) becomes ^-independent and the linear in 
B term in Eq. (|71l) is absent. Thus, in this case the su- 
perconductivity can be promoted by a parallel magnetic 
field, at least in the weak field limit. This agrees with the 
results of Ref . j3j| , where the gradient term linear in B 
and D was introduced on the phenomenological grounds 
for a surface superconductor. The order parameter which 
occurs at T c at non-zero B is modulated in the xy plane: 
, 7 (r) = r] e iQr , with Q (x (z x B) [13, see also Ref. [24| . 
It should be noted though that the field-induced increase 
in T c may indicate the onset of a magnetic instability of 
the superconducting state, the investigation of which is 
beyond the scope of the present work. 



V. CONCLUSIONS 

We studied the magnetic properties of a clean super- 
conductor with spin-orbit coupling. We focussed on the 
weak-field limit near the critical temperature, where the 
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APPENDIX A: DERIVATION OF EQ. (IMl) 

To derive the free energy for a nonuniform distribution 
of the order parameter, we start with a representation of 
the partition function for the BCS Hamiltonian (|27l) in 
terms of a functional integral over the Grassmann fields 
Cfe Q (r) and c fcQ (r): 



Z = / VcVc e' 



(Al) 



where S = J Q dr [J2k c kad T Cka + H(t)]. The interaction 
term in the action can be written as 

Smt = ~E f drJ2Bt(q,r)B a (q,T), 



where 



-fe+q/2,a ( T )Cfc+q/2,/3 ( T )- 

k 

The interaction term is then decoupled by means of 
the Habbard-Stratonovich transformation, introducing a 
complex bosonic field r/ a (q, r): 

e~ s ^ - J V V * a V Va exp{-^ jy^^al 2 
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1 [ p 

/ dT ^2 A "/3( fc '9;T)Cfc+q/2 : a(T)c_ fc+ q/2, /3 (T) 
fcq 



The last two terms in the exponent can be written as 

kq 

+H.c, 

where 

A Q(3 (fc,q;r) = ^ ^(9, r)* a>a/3 (fe) (A2) 
a 

is the order parameter matrix in the pseudospin space 
[cf. Eq. (EJl]. 

The next step is to integrate out the fermionic 
degrees of freedom, which can be achieved by us- 
ing the four-component Nambu spinor fields Cfc(r) = 
[cfca(r), C-k a (T)] T and calculating a Gaussian fermionic 
integral. As a result we arrive at the following represen- 
tation of the partition function: 



Z = J V V * a V Va e- s 'ff ["*>"], 



(A3) 



where 



S *lf = ^Yl J dr^2\ Va \ 2 -^Trln(l-go^) (A4) 



with the order parameter matrix defined by Eq. JA2). 
The trace in the action l|A4Jl should be understood as 
the matrix trace in the four-dimensional Nambu x pseu- 
dospin space, accompanied by the operator trace in the 
fcr-space. 

Using the partition function (jA.Sfl , we can calculate the 
free energy of the system: T — —(l/(3)\nZ. The BCS 
mean-field approximation corresponds to a stationary 
saddle point of the effective action (|A4|) . For rj a (q,T) — 
r]a(q), the saddle-point action becomes 5^// = f3!F, with 
the free energy (or, more precisely, the difference between 
the free energies of the superconducting and the normal 
states at the same temperature) given by 



T = v E E i^' 2 - M Trln(1 - 0oE) - (A7) 

a q 

The order parameter components satisfy the saddle-point 
equations 5 J- /Srj* = (the GL equations). In the vicinity 
of the critical temperature at arbitrary magnetic field, 
the order parameter is small, so we can keep only the 
quadratic in r/ a terms in the expansion of the trace in the 
free energy l|A7|l . In terms of the Fourier-transformed 
basis functions 



is the effective action for the superconducting order pa- 
rameter. Here Go is the Gor'kov-Nambu Green's function 
at rj = rf — (i.e. in the normal state): 



*a,a/3(p) 



E 

k 



(A8) 



Go 



G 

-G T 



(A5) and the Green's functions JSJ), we have 



where G = (—d T — £) _1 is a 2 x 2 matrix in the pseu- 
dospin space, which satisfies Eq. JJJ), and £ is the 4x4 
matrix self-energy function describing the superconduct- 
ing pairing: 



F = Y1 J '^idr 2 <(ri)5 o6 (ri,r 2 )77 6 (r 2 ), (A9) 

ab 



A 

At 



(A6) 



with the kernel 



S a b{ri, r 2 ) = y5ij5(ri - r 2 ) 

-\ T H I K*pfa) G fh (*i + Y' r 2 + Y ;uJn ) *wteO G «rf ( ri - f' r2 - Y ; ^ n ) • (A10) 



Substitution of the factorized Green's function I|1U|) in (|A10ll gives the phase factor 



exp 



lip 



( , Pi . P2\ . . / Pi P2 



exp 



2i(p{r u r 2 ) 
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[to prove this, one can use the Taylor expansions of the 
yj's with respect to P1.2, and also the identities (| 1 1|) ] . The 
next step is to use 

cxp ji^ J 2 A(r)dr| r/(r 2 ) = e - i ( ri - P2 > Dl 77(ri), 



where = — iV r + (2e/Hc)A, to cast the free energy 
(|A9() in the form with the function S(R) given by 



S a6 (il) = \ T ^2 j d Pl d P2 *l,a^(Pl)*b,75(p2)exp 

Pi - P2 



xG 



01 



R 



r 



*^B(P1 x p 2 ) 



Pi - P2 



(All) 



Finally, taking the Fourier transform of this expression, 
we arrive at Eq. 

The analysis in the non-centrosymmetric case can be 
done in a similar fashion, the only difference being that 
there is no pseudospin degrees of freedom, and G, 'J', and 
A become just scalar functions. The partition function 
still has the form i|A3|l , but the effective action now reads 



S eff = ^J2 1 ^^|, ?a | 2 -iTrln(l-g S), (A12) 



where Go an d £ are 2x2 matrix operators in the Nambu 
space and the fcr-space. Repeating all the steps leading 
to Eq. (|A1 If) , we arrive at Eqs. <|55jl and 



APPENDIX B: GRADIENT ENERGY NEAR T c 

In this Appendix we estimate the lowest eigenvalue 
of the matrix differential operator JC a t = K a bijDiDj, 
where K a \,^ are constant coefficients, a, b = l..d, and 
i, j = x, y, z. We choose B along the z axis, i.e. B = Bz 
(one can always achieve that by rotating the coordinate 
system, which is equivalent to a re-definition of A" a {,.y ). 
It is convenient to introduce new operators 



(Bl) 



It is easy to check that the operators a± satisfy the rela- 
tions a + — a! and [a_,a + ] = 1, and therefore have the 
meaning of the lowering and the raising operators respec- 
tively, while the operator 03 = 0,3 commutes with both 
of them: [03, a±] = 0. Representing YL a b in terms of the 
operators (|B1|) . we have 




fCab 



eB_ 

he 



where the coefficients K a b t . nm are linear combinations of 
K a b,ij and therefore do not depend on B. It immediately 
follows from the last expression that all eigenvalues of JC 
are linear in B. 

To calculate the eigenvalues explicitly, it is convenient 
to choose the basis of states \N, p) such that 



a + \N,p) = VWTT\N +l,p) 
a-\N,p) = Vn\N-1,p) 
a 3 \N,p)=p\N,p), 

where N = 0, 1, ... has the meaning of the Landau level 
index and p is a real number which is proportional to 
the wave vector along the z-axis: p — k z yJhc/eB. Ex- 
panding the eigenfunctions of JC in this basis: rj a {r) — 
J2n p C a ,N ,p{ r \N , p) , we arrive at a system of linear equa- 
tions for the coefficients C a _N,p, which is infinite in gen- 
eral. The upper critical field then corresponds to the min- 
imum eigenvalue of this system with respect to p (while 
it is usually assumed that the minimum is achieved for 
p = 0, some exceptions are discussed, e.g. in Ref. 0). 

In some simple cases, the diagonalization procedure 
outlined above can be carried out analytically. For exam- 
ple, for a one-component order parameter in an isotropic 
s-wave superconductor we have 



t = K{Dl + Dl + Dl) = ^- 



-K(4a + a-+al + 2). (B3) 
Since a + a_ \N,p) — N\N,p), we have 



E 

n.m — ±,3 



t\N,p) = ^-K{AN + p 2 + 2)\N,p). (B4) 
nc 

The lowest eigenvalue corresponds to N = p = 0, which 
gives the standard expression for the critical temperature 
suppressed by the field: 

T c (B)=T c (B = 0)-^-B. (B5) 
he a 



K-ab,nmQ> n @"m •> 



(B2) 
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